CRREL 


For  conversion  of  SI  metric  units  to  U.S. /British 
customary  units  of  measurement  consult  ASJM 
Standard  E380,  Metric  Practice  Guide,  published 
by  the  American  Society  for  Testing  and  Materials, 
1916  Race  St.,  Philadelphia,  Pa.  19103. 


Cover:  Large  snowdrift  on  the  roof  of  a school  in 

Buffalo,  New  York.  (Photograph  by  Robert 
Red  field.) 


CRREL  Report  77-12 


Roof  loads  resulting  from  rain-on-snow 

Samuel  C.  Colbeck 


j May  1977 


Unclassified 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  D»t»  Entered) 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


REPORT  DOCUMENTATION  PAGE 


3.  RECIPIENT'S  CATALOG  NUMBER 


12.  GOVT  ACCESSION  NO. 


CRRELft*** -77-12 


5 TYPE  OF  REPORT  A PERIOD  COVERED 


ROOF  LOADS  RESULTING  FROM  RAIN-ON-SNOW. 


6 PERFORMING  ORG.  REPORT  NUMBER 


AUTHORfi; 


Samuel  C/fcolbeck 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

U.S.  Army  Cold  Regions  Research  and  Engineering  Laboratory 
Hanover,  New  Flampshire  03755 


DA  Projcct/4A16l  101A91D 
Task  256  * 


It.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 


U.S.  Army  Cold  Regions  Research  and  Engineering  Laboratory 
Hanover,  New  Hampshire  03755 

4.  MONITORING  AGENCY  NAME  A AOORESSfff  different  from  Controlling  Office) 


15.  SECURITY  CLASS,  (of  this  report) 


Unclassified 


1 5a.  DECLASSI  F| CATION/ DOWNGRADING 
SCHEDULE 


16.  DISTRIBUTION  STATEMENT  (of  this  Report) 


Approved  for  public  release;  distribution  unlimited 


17.  OlSTR\BUT\ON  STATEMENT  (of  tha  abstract  entered  in  Block  20,  different  from  Report) 


18.  SUPPLEMENTARY  NOTES 


19.  KEY  WORDS  ( Continue  on  reverse  aide  if  necessary  and  identify  by  block  number) 

Drainage 
Loads  (forces) 

Rainfall 

Roofs 


ABSTRACT  (Continue  on  reverse  side  If  necessary  and  Identify  by  block  number) 

A computer  program  to  calculate  the  increased  live  load  on  a snow-covered  roof  due  to  rain-on-snow 
the  25-year  rainstorm  (ailing  on  a heavy  snow  load  on  a flat  roof  in  Hanover,  New  Hampshire,  an  adc 
(20  lb/ft2)  of  liquid  water  is  added  to  the  live  load.  The  additional  load  due  to  rain-on-snow  is  very 
snow  properties  and  characteristics  of  the  roof.  A wide  range  of  live  loads  is  possible,  depending  on 
circumstances. 


FORM 
I IAN  7J 


EOITION  OF  t NOV  »S  IS  OBSOLFTE 


SECURITY  CL  ASSfFlC  ATlOw  O F THIS  PAGE  ''Wtfn  Data  Entered: 


PREFACE 


This  report  was  prepared  by  Dr.  Samuel  C.  Colbeck,  Geophysicist,  of  the  Snow 
and  Ice  Branch,  Research  Division,  U.S.  Army  Cold  Regions  Research  and  Engi- 
neering Laboratory.  The  research  was  funded  by  DA  Project  4AI 61  101  A91  D, 
In-House  Laboratory  Independent  Research,  Task  256,  Roof  Loads  Resulting 
front  Rain-on-Snow. 

Dr.  Haldor  Aamot,  Thaddeus  Johnson  and  Robert  Redtield  of  CRREL  techni- 
cally reviewed  the  manuscript.  Their  contributions  to  this  project  aie  gratefully 
acknowledged. 

I • 


S 


CONTENTS 


i 


i 

> 


t 


; 


\ 


I 


K 


w 

f 

» 


.M 


Page 


Abstract i 

Preface ji 

Summary v 

Nomenclature vi 

Introduction I 

Water  movement  through  snow  1 

Vertical  percolation  1 

Lateral  flow  2 

Plat  roofs 3 

Sloping  roofs 3 

Rainfall  intensity-duration  effects 7 

Miscellaneous  effects  9 

Radial  flow  to  drains 9 

Plow  along  gutters  on  snow-covered  roofs  1 3 

Snow  structure 14 

Basal  layer 14 

Conclusion 14 

Literature  cited 15 

Appendix  A.  Computer  program  calculating  roof  loads  from  rain-on-snow  17 


ILLUSTRATIONS 


I igure 

1 . Snow-covered  flat  roof  with  an  internal  drain I 

2.  Roof  load  due  to  refro?en  rainwater  shown  as  a function  of  snow 

depth  for  various  snow  temperatures 2 

3.  Thickness  of  the  saturated  layer  at  the  root  margin  shown  as  a func- 

tion of  time  for  the  specified  conditions 4 

4.  Roof  with  a gradual  slope  4 

5.  Dimensionless  volume  of  water  in  the  saturated  layer  shown  for 

given  parameters  as  a lunction  of  roof  slope 5 

6.  Inverse  of  the  correction  factor  F shown  as  a function  of  slope 

for  n = 1 0 6 

7.  Slope  correction  factor  shown  as  a function  of  slope  for  various 

values  of  ( ak/l),/'‘  6 

8.  Volume  of  water  per  unit  width  of  saturated  layer  shown  as  a 

function  of  time  for  two  roof  slopes 7 

9.  Total  live  load,  fraction  ol  rainfall  retained,  and  depth  of  water  im- 

pounded on  a flat,  two-dimensional  roof  for  the  25-year  rainfall 
at  Hanover,  N. 1 1 8 


iii 


I 


k 


Figure 

10.  Saturated  layer  of  snow  on  a flat  roof  shown  for  the  cases  where 

the  saturated  layer  does  and  does  not  cover  the  drain 

1 1 . Three  advantages  of  a recessed  drain  over  a flush  drain  

I 2.  Profile  of  the  saturated  layer  shown  as  a function  of  R for  various 

drain  sizes  for  radial  flow  over  a flat  roof 

13.  Ratio  of  the  volume  of  water  in  the  saturated  layers  of  radial  and 

two-dimensional  roofs  shown  as  a function  of  the  drain  size  on  a 
radial  roof  

14.  Flow  on  a flat  roof  for  two  conditions 

15.  Depth  of  the  saturated  layer  shown  as  a function  of  time  for  an  in- 

tense rainstorm  for  two  values  of  permeability  of  the  basal  layer. 


SUMMARY 


Rain  tailing  on  a heavy  snow  load  can  cause  roof  failures  by  adding  temporary 
weight  to  a roof.  This  additional  weight  can  be  calculated  as  a function  of  time 
for  any  given  depth,  porosity,  temperature  and  permeability  of  snow;  size,  slope, 
and  shape  of  roof;  and  duration  and  intensity  of  rainstorm.  When  the  25-year 
rainstorm  for  Hanover,  New  Hampshire,  is  used  as  an  example,  it  is  found  that  a 
two-dimensional,  flat  root  could  retain  enough  rainwater  to  reach  about  50%  of 
most  design  live  loads.  This  weight  could  be  partially  reduced  by  giving  the  roof 
a slight  inclination;  however,  the  weight  would  be  increased  if  the  water  moved 
radially  to  internal  drains  rather  than  to  large  gutters. 

It  is  shown  that  the  roof  load  due  to  rain-on-snow  is  very  sensitive  to  spacing 
between  drains,  snow  depth  and  roof  inclination.  The  load  on  a flat  roof  is  less 
sensitive  to  the  duration  of  the  25-year  rainstorm  and  the  permeability  of  the 
snow.  Roofs  with  radial  flow  to  a single,  central  drain  tend  to  drain  more  slowly 
than  roofs  with  parallel  flow  to  large  gutters.  The  advantages  of  enlarged,  re- 
cessed drains  are  shown. 

In  the  Appendix,  a computer  program  is  given  for  calculating  the  total  weight 
of  wet  snow  as  a function  of  duration  for  any  design  basis  rainstorm.  In  addition 
to  choosing  the  design  basis  rainstorm,  the  snow  depth,  porosity,  temperature, 
and  permeability,  as  well  as  the  roof  size,  slope,  and  shape  must  be  specified  by 
the  user. 
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d depth  of  saturated  layer 
dL  depth  of  saturated  layer  at  drain 
dQ  maximum  depth  of  saturated  layer 
D grain  size 
F slope  correction  factor 
g gravitational  constant 
h thickness  of  snow 
/ rainfall  intensity 
/ infiltration  into  the  saturated  layer 
&s  intrinsic  permeability  of  saturated  layer 
ku  intrinsic  permeability  of  unsaturated  layer 
L roof  length 
Lc  culvert  length 
n \fakjT 

q discharge  through  the  saturated  layer 
C?D  drain  discharge 
r radial  coordinate 
/•R  drain  radius 
re  roof  radius 
R rjre 

rotre 
R,  X/L 

5W  water  saturation  as  percent  of  pore  volume 
Swi  irreducible  water  saturation 
5*  (Sw-Swi)/(1 -5wi) 

time 

T snow  temperature  (°C) 

u volume  flux,  water  volume  per  unit  area  per  unit  time 
V water  volume  of  saturated  layer 
w culvert  width 
W total  weight  per  unit  area 
Wu  liquid  weight  per  unit  area  in  the  saturated  layer 
x coordinate  direction  along  roof 

, speed  of  a value  of  S* 
t/r  | S*  p 

a 5.47  x 10''  nr1  s’1 

0 roof  slope 
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ROOF  LOADS  RESULTING  FROM  RAIN-ON-SNOW 


Samuel  C.  Colbeck 


INTRODUCTION 


WATER  MOVEMENT  THROUGH  SNOW 


The  cost  of  constructing  a modern  building  is  so 
great  that  much  time  and  effort  must  be  devoted  to 
saving  material  and  labor.  Unfortunately  the  efficiency 
of  the  design  is  often  limited  by  the  many  variables 
which  must  be  considered.  These  variables  including 
winds,  soil  conditions,  and  snow  loads  — are  often  dif- 
ficult to  quantify  and  therefore  difficult  to  accommo- 
date in  the  design.  One  important  aspect  of  many 
designs  roof  load  due  to  rain  falling  on  a snow- 
covered  roof  - is  described  here. 

Much  consideration  has  been  given  to  regional  snow 
accumulations,  including  both  depths  and  densities 
(e.g.  Tobiasson  and  Redfield  1973).  Thorburn  and 
Schriever  (1959)  discussed  roof  failures  resulting  from 
rain-on-snow,  emphasizing  the  significance  of  rain 
falling  on  an  already  heavy  snow  load.  Rain  has  two 
special  properties  in  this  regard.  First,  rain  is  much 
more  dense  than  snow  but  can  be  added  to  the  snow 
without  increasing  its  bulk.  In  fact,  an  observer  may 
believe  incorrectly  that  the  roof  load  decreases  as  the 
snow  densifies  following  the  introduction  of  rain  water. 
Second,  rain  falling  on  snow  percolates  downward  until 
it  reaches  the  roof  and  then  spreads  laterally  or  flows 
downslope.  If  the  snow  is  at  subfreezing  temperatures, 
some  of  the  percolating  water  will  be  refrozen  in  place, 
but  if  the  snow  is  already  wetted,  all  the  water  will  pass 
through  the  snow  and  drain  off  the  roof.  In  the  former 
case  some  of  the  water  will  remain  on  the  roof  as  an 
integral  part  of  the  snow  load,  whereas  in  the  latter  case 
all  the  water  passes  through  the  snow,  adding  only 
transient  weight  to  the  roof.  Our  purpose  here  is  to 
establish  techniques  for  calculating  this  weight  as  a 
function  of  time  for  any  given  set  of  conditions.  The 
basic  question  is  how  much  of  the  transient  rainwater 
can  be  retained  by  the  snow  on  the  roof  at  any  time. 

For  a "design  basis  rainstorm,"  the  roof  load  is  con- 
sidered to  be  at  a maximum  and  it  is  this  load  which 
the  roof  must  be  designed  to  support. 


Rainwater  falling  on  a snow-covered  roof  percolates 
downward  until  it  reaches  the  impermeable  roof  surface. 
The  water  then  flows  along  the  roof  surface  to  an  inter- 
nal drain  or  the  roof  edge.  The  two  modes  of  flow 
unsaturated  vertical  percolation  and  saturated  lateral 
flow  over  the  roof  - are  quite  distinctive  and  require 
separate  descriptions  (see  Fig.  1). 
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Figure  1.  Snow-covered  flat  roof 
with  an  internal  drain  where  h is 
thickness  of  snow,  L roof  length, 
d0  maximum  depth  of  saturated 
layer,  and  x coordinate  direction 
along  the  roof.  The  saturated 
layer  of  water  flowing  over  the 
roof  is  elliptical  in  shape. 


Vertical  percolation 

The  condition  of  snow  prior  to  the  rainfall  is  very 
important  because  of  the  wide  variety  of  responses  that 
can  happen  in  any  given  situation.  If  the  snow  is  at  a 
subfreezing  temperature,  some  of  the  rainwater  will  be 
frozen  in  place  in  order  to  raise  the  snow  to  the  melting 
point.  This  volume  per  unit  area  is  equal  Xo  - p^Th/]60 
where  ps,  T,  and  h arc  the  density,  temperature  and 
depth  of  the  snow,  respectively.  As  shown  on  I igure  2 
for  a typical  snow  density  of  0.5  Mg/m'1 , the  thermal 


Figure  2.  Root  loud  due  to  /e  frozen  rainwater 
shown  as  a function  of  snow  depth  for  various 
snow  temperatures.  The  snow  density  is  assumed 
to  be  0.5  Mg/m3. 

el  feci  adds  only  3. 1 25  kg/m2  of  load  (0.64  I b/ 1 1 2 ) 
per  meter  of  depth  and  degree  Celsius  of  subfreezing 
temperature.  It  seems  unlikely  that  this  effect  could 
be  more  than  a small  percentage  of  the  total  roof 
load  because  heavy  rains  of  sufficient  duration  are 
very  unlikely  when  a thick,  cold  snowload  occurs. 

The  unsaturated  percolation  of  water  through  the 
upper  portion  of  the  snow  also  adds  to  the  roof  load 
because  of  the  need  to  increase  the  liquid  water  con- 
tent in  order  to  allow  water  percolation.  If  the  flow 
through  the  snow  is  uniform  and  homogeneous  (i.e. 
does  not  flow  in  distinct  drainage  channels),  this  in- 
crease in  weight  due  to  the  liquid  saturation  can  be 
calculated  from  known  principles  of  the  unsaturated 
water  flow  in  snow  (e.g.  Colbcck  and  Davidson  1973, 
Colbeck  1976).  The  intensity  of  the  rainfall  and  the 
grain  size,  density  and  depth  of  the  snow  would  have 
to  be  known  to  make  this  calculation.  Assuming  a 
constant  source  of  rain  of  sufficient  duration  to 
penetrate  the  entire  snowcuver,  the  weight  of  liquid 
per  unit  area  would  be  equal  to 

=/tpw0[(/Mu)'/3(1  -Swi)+Swi|  (1) 

where  pw  is  density  of  water,  0 snow  porosity,  / in- 
tensity of  the  rain,  a a constant  (5.47  x 106  rtf'  s'1 ), 
and  Swi  the  irreducible  water  saturation  (about  0.07 
depending  on  the  degree  of  snow  metamorphism).  A’u 
is  the  intrinsic  permeability  given  by  Shimizu  (1970)  as 

hu  =0.077  D2  exp  (-7.8  ps)  (2) 

where  D is  the  grain  size.  For  a rainstorm  of  10'6 
m3/m2  s (~  0.14  in./h)  with  a snow  depth  of  1 m and 


permeability  of  tO'10  m2,  the  inctcasc  in  tool  load  due 
to  the  moving  water  would  be  76  kg/m2  (15.5  lb/ft2) 
per  meter  of  depth.  This  weight  alone  approaches  half 
of  most  live  load  design  limits  if  the  rainstorm  lasts  long 
enough  to  deposit  this  much  rain  and  if  the  llow  is  uni 
form  throughout  the  snow.  In  fact  neither  of  these 
conditions  is  likely.  As  shown  later,  the  probability  of 
a rainstor m of  this  intensity  and  duration  falling  on  snow 
is  small,  furthermore,  the  roof  load  would  probably  be 
reduced  somewhat  below  the  value  calculated  here  be- 
cause of  nonuniform  flow  through  the  snow.  I he  de- 
gree of  nonuniformity  would  be  highly  variable,  de- 
pending on  the  nature  ol  the  snow  on  the  roof  prior  to 
the  onset  of  the  rain.  Thus,  while  some  correction  could 
be  applied  to  account  for  the  channeling  often  observed 
in  snow  on  the  ground,  there  is  no  way  of  assessing  what 
this  correction  factor  might  be.  Since  it  is  conservative 
to  assume  a uniform  distribution  of  the  flow,  the  cal- 
culations are  made  accordingly. 

The  effects  of  the  intensity  of  the  rainfall  and  perme- 
ability of  the  snow  are  reduced  by  their  one-third  power 
dependence.  An  order-of-magnilude  increase  in  rainfall 
intensity  only  about  doubles  the  weight  due  to  the 
unsaturated  snow.  As  shown  later  the  more  intense 
rainfalls  are  generally  shorter  in  duration;  thus,  the  less 
intense  but  longer  lasting  rainfalls  can  add  more  weight 
to  the  roof  than  the  more  spectacular  (but  shorter) 
heavy  rains.  Furthermore,  the  less  intense  rainstorms 
are  more  likely  to  cause  a homogeneous  soaking  of  the 
snow;  hence,  rapid  runoff  due  to  the  formation  of  flow 
channels  is  less  likely  for  rainstorms  of  lower  intensity . 

Grain  growth,  density  increase,  and  a general  homog- 
enization of  the  snowcover  occur  when  large  quantities 
of  rain  or  meltwater  first  infiltrate  the  snowcover.  The 
properties  of  the  snowcover  change  rapidly  during  this 
period  of  “melt  metamorphism”  and  the  permeability 
can  increase  significantly.  Again  because  of  the  one- 
third  power  dependence  of  weight  on  permeability,  the 
roof  load  due  to  the  unsaturated  snow  would  decrease 
very  slowly  with  increasing  permeability.  Because  of 
the  general  lack  of  observations  of  the  physical  proper- 
ties of  snow  on  roofs,  it  is  difficult  to  assign  values  to 
the  intrinsic  permeability  of  the  snow  but  the  assumed 
value  of  10'10  m2  is  in  the  correct  range  for  the  un- 
saturated snow. 

Lateral  flow 

Once  the  percolating  water  reaches  the  roof  surface, 
the  water  moves  laterally  in  a saturated  layer.  The 
permeability  to  the  water  in  this  saturated  layer  is 
orders  of  magnitude  greater  than  in  the  unsaturated 
layer  for  two  reasons.  First,  the  permeability  to  the 
water  increases  as  the  cube  of  the  liquid  saturation 
(Colbcck  and  Davidson  1973)  and  the  unsaturated 


(6) 


layer  is  ty  pically  1 0,  saturated  while  the  saturated 
layer  is  nearly  100%  saturated.  Thus  the  saturation 
effect  increases  the  permeability  to  the  liquid  in  the 
saturated  layer  by  a factor  ol  about  103.  Second, 
the  grain  si/e  in  the  saturated  layer  is  likely  to  he 
greater  because  ol  “'temperature  gradient  meta- 
morphism” in  the  warmer  basal  layer  prior  to  the 
introduction  ol  liquid  water  and  because  ot  the  rapid 
grain  growth  in  saturated  snow  alter  the  introduction 
of  liquid  water  (Wakahama  1968).  As  shown  by  eq  2, 
the  intrinsic  permeability  increases  as  the  square  of 
the  grain  size.  Accordingly,  an  intrinsic  permeability 
of  lO"1'  mr  might  be  typical  tor  the  saturated  layer 
A-s. 

At  the  root  surface  the  pressure  gradient  available 
to  drive  the  lateral  flow  is  very  small  and  the  flow  is 
reduced  accordingly,  f or  a house-size  gable  roof 
of  10°  slope,  the  weight  of  the  saturated  layer  could 
reach  about  !0  kg.'mr  (2  Ib/ft 2 ) assuming  no  blockage 
at  the  rool  edge.  Since  most  gable  roofs  are  steeper 
than  this,  the  weight  of  their  saturated  layers  would 
generally  be  negligible.  At  slopes  of  less  than  10u,  the 
flow  is  reduced  significantly  and  the  thickness  and 
weight  of  the  saturated  layer  are  increased  accordingly. 
Lateral  flow  on  both  flat  and  low-slope  roofs  is  de- 
scribed below. 

Flat  roofs 

I or  the  idealized,  two-dimensional  roof  shown  in 
Figure  I , the  profile  of  the  saturated  layer  is  often 
assumed  to  be  elliptical,  l or  steady  flow  where  no 
blockage  occurs  at  the  drain,  the  profile  is  given  by 

(£\2  + (*-—*11  =1.  (3) 


where  d is  the  depth  of  the  saturated  layer. 

The  volume  of  water  V stored  in  the  saturated  layer 
per  unit  width  of  roof  is 

V - 0.5  nL~  s/l/ak^Q  (4) 

where  / is  the  flow  into  the  saturated  layer  from  the 
overly  ing  unsaturated  layer  (m  ’ m2  s),  0 is  snow 
porosity  (I  - pjp j),  and  is  intrinsic  permeability  of 
the  saturated  layer.  1 he  maximum  thickness  of  the 
saturated  layer  (t/0)  for  steady  flow  is 


This  equation  is  solved  here  by  assuming  that  transient 
wave  elfects  are  small  and  that  the  elliptical  profile  is 
always  maintained.  I hen  the  equation  can  be  simplified 


This  equation  can  be  immediately  solved  when  1(c)  is 
either  a constant  or  zero,  c.g.  during  the  onset  and  then 
abrupt  cessation  of  rain.  The  solutions  for  these  cases 
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where  c/0  (0)  is  the  value  of  d0  when  / changes.  An  ex- 
ample of  dQ  (l)  is  given  on  Figure  3.  The  thickness  of 
the  saturated  layer  increases  rapidly  following  the  onset 
of  infiltration,  although  it  is  very  unlikely  that  more 
than  50%  of  the  maximum  possible  thickness  would 
ever  be  realized  because  such  heavy  rains  are  not  of 
sufficient  duration.  Once  the  infiltration  ceases,  the 
depth  of  the  saturated  layer  decreases  rapidly  until 
more  rain  falls  on  the  snow. 

The  total  snow  load  per  unit  width  carried  by  a flat 
roof  with  an  internal  drain  is  equal  to  the  weight  of  the 
snow  plus  the  liquid  in  the  unsaturated  layer  plus 
O.5rr0/.2pw  (//q*J|/j,  the  weight  of  the  liquid  in  the 
saturated  layer  from  eq  4.  The  saturated  layer  is  some- 
what more  sensitive  to  changes  in  infiltration  or 
permeability  than  the  unsaturated  layer  although,  be- 
cause of  the  square  root  dependence,  the  weight  of  the 
saturated  layer  is  less  responsive  to  As  than  might  have 
been  thought.  At  small  values  of  time,  the  As  depend- 
ence of  the  saturated  layer  nearly  vanishes  as  shown 
below.  Note  that  flat  roofs  are  very  sensitive  to  the 
sparing  of  the  drains  since  the  weight  per  unit  width  of 
saturated  layer  increases  as  L* . 

Sloping  roofs 

The  complete  equation  describing  flow  along  a 
sloping  roof  (Colhcck  1974)  is 


l - e nonsteadv  flow  where  / varies  with  time,  the  fol 

lowing  equation  (Colbcck  1973)  describes  the  deptfi 

of  the  saturated  layer  as  functions  of  v and  I: 


*"!r  " 


where  6 is  the  roof  slope. 


* - ■ ' 


Time  (s) 


Figure  3.  I hickness  of  the  saturated  layer  at  the  roof  margin  shown  as  a func- 
tion of  time  for  the  specified  conditions.  The  thickness  increases  rapidly  fol- 
lowing the  onset  of  rain  but  it  is  unlikely  that  such  a heavy  rain  could  last  long 
enough  to  generate  more  than  50%  of  the  maximum  possible  value  of  dQ. 

Once  the  infiltration  due  to  the  rainfall  stops,  the  thickness  of  the  saturated 
layer  decreases  rapidly. 


I 


Figure  4.  Roof  with  a gradual  slope.  For  a slightly  dip- 
ping roof  the  thickness  of  the  saturated  layer  is  inter- 
mediate between  a flat  and  a steeply  dipping  roof. 

For  a steeply  dipping  roof  this  equation  can  be  simpli-  Hence  for  steady  flow, 

fied  by  neglecting  the  second  order  derivative,  whereas 

for  a flat  roof,  the  first  term  is  zero.  Roofs  with  /.  a aW  \ , , ,,,, 

gradual  slopes  are  also  common,  and  il  is  very  impor-  \ dx  / 

tant  to  analyze  them  as  accurately  as  possible.  Because 

of  the  extreme  sensitivity  of  flow  to  slopes  at  angles  where  0 is  about  equal  to  sin  0 for  shallow  slopes.  Al- 

of  less  than  10°,  either  of  the  simplifications  mentioned  though  eq  1 1 has  no  temporal  dependence,  it  cannot 

above  could  cause  large  errors.  As  for  flat  roofs,  we  be  readily  solved  because  of  its  nonlinearity.  Fortu- 

begin  by  establishing  the  upper  limit  of  weight  on  a natcly  we  are  not  interested  in  solving  lor  d per  se  but 

roof  of  gradual  slope  (see  Fig.  4)  during  a steady  rain.  rather  we  can  solve  for  the  volume  of  water  1/  in  the 

The  flux  u (volume  per  unit  area  per  unit  time)  along  saturated  layer: 

the  roof  is  given  by 

L 

u = a*s  (sin  0 - j . (10)  V = <t>J  d dx  . (12) 

0 
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- Exact  solution  for  flat  roof 
for  v'ohTI  - 10 


Approximate  Solution 

for  v ak/  I =10 


Equation  (14) 
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Figure  5.  Dimensionless  volume  of  w a ter  in  the  saturated  layer  [(2otkJ 
0//- 2)  V / shown  for  given  parameters  as  a fun c tion  of  root  slope.  Die 
solution  for  V given  by  eq  14  is  only  approximate  and  as  the  root  slope 
decreases,  the  error  increases.  Accordingly  an  approximation  to  the 
solution  is  given  tor  small  values  of  0. 


From  eq  I ! and  I 2 we  get 


I (dl-dl  + la'  k;'  L2)  (13) 


where  dL  is  the  depth  of  the  saturated  layer  at  the 
drain.  Lquation  1 3 is  an  exact  expression  for  the 
volume  of  water  in  the  saturated  layer,  but  the  bounda 
ry  values  J|  and  d0  cannot  be  specified  a priori.  This 
equation  could  be  used  with  (dj  -</q)  as  a parameter, 
but  for  the  practical  purposes  ol  deciding  how  much 
water  can  be  impounded  by  the  saturated  layer,  we 
must  make  the  assumption  that 

/«•'*;'  ‘-2  K-^|- 

Accordingly , eq  13  can  be  approximated  by 


slope  on  Figure  5.  There  is  a small  error  in  the  solution 
for  slightly  sloping  roofs  because  eq  14  tends  to  over- 
estimate V at  small  slopes.  This  error  arises  because  we 
have  neglected  d £ - df,,  a negative  number.  An  approxi- 
mate solution  for  the  volume  of  water  in  the  saturated 
lay  er  on  a shallow  root  is  also  shown  on  Figure  3.  The 
approximation  is  made  by  noting  that  for  a value  ol 
(a k I)  equal  to  10,  l’  does  not  apnroach  inf  inity  as 
0 vanishes  as  suggested  by  eq  14,  but  rather  T has  an 
upper  limit  given  by  eq  4 tor  flat  roots.  Also,  as  the 
slope  increases,  the  actual  solution  must  approach  the 
approximate  solution  of  eq  14.  The  difference  between 
eq  14  and  the  approximate  solution  is  only  significant 
at  angles  below  10  where  the  root  load  is  most  sensi- 
tive to  slope. 

To  include  the  ellect  ol  roof  slope,  we  deline  a slope 
correction  factor  / as 


V ■ !,-  00-'  la  ' I2 


which  is  the  same  as  the  formula  lor  the  volume  ol 
water  on  a steeply  sloping  roof,  because  the  assump- 
tion stated  above  is  only  valid  lor  sleeper  slopes. 

The  volume  ol  water  contained  in  the  saturated 
layer  as  given  by  eq  1 4 is  shown  as  a function  of  root 


I (())  (1 

' »n(0)  11 

where  K'n(0)  is  the  weight  ol  the  saturated  laser  for 
any  slope  I)  or  parameter  »,  where 

n \fakji . (I 
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Figure  S.  Volume  of  water  per  unit  width  of  satu- 
rated layer  shown  as  a function  of  time  for  two  roof 
slopes.  I he  difference  between  the  two  slopes  increases 
with  time  until  40  x I O2  s when  the  percolating  water 
stops  reaching  the  roof.  The  more  gradual  roof  carries 
about  twice  the  weight  of  the  steeper  roof  at  40x1  O ' s. 

Figure  7 shows  clearly  that  the  effect  of  increasing 
the  roof  slope  in  order  to  provide  better  drainage  is 
highly  dependent  on  the  values  of  (akJI)Vl . The  weight 
of  the  ponded  layer  of  water  decreases  sharply  for  small 
increases  in  slope  at  slopes  of  less  than  10°,  but  the  opti- 
mum slope  is  dependent  upon  the  values  of  rainfall 
intensity  and  snow  permeability  chosen  for  the  optimi- 
zation calculation,  For  both  gently  sloping  and  flat 
roots,  the  weight  of  the  saturated  layer  per  unit  width 
increases  as  the  value  of  the  spacing  between  the  drains 
squared  4 L- . 1 his  is  a significant  result  which  shows 
the  critical  value  of  drain  spacing  on  shallow  roofs. 

Figure  8 shows  die  volume  of  water  per  unit  width 
of  the  saturated  layer  as  a function  of  time  for  two 
slopes  where  infiltration  stops  after  40x103  s.  The 
more  shallow  slope  carries  a larger  volume  of  water 
but  it  takes  more  time  to  reach  its  maximum  value. 
During  the  first  hour  of  infiltration  into  the  saturated 
layer,  the  ef  fect  of  slope  is  not  as  significant  as  sug- 
gested by  Figure  5.  The  effect  of  slope  increases  as 
infiltration  continues,  until  ultimately  the  more  shallow 
slope  carries  more  than  twice  as  much  water  in  its 
saturated  layer  as  does  the  steeper  slope.  Also,  the 
more  shallow  slope  responds  more  slowly  once  the  infil- 
tration ceases  by  taking  longer  to  drain . Flowcver,  be- 
cause it  is  very  unlikely  that  such  an  intense  rainstorm 
will  last  for  this  period  of  time,  the  difference  due  to 
roof  slope  is  reduced  somewhat.  The  decreased  response 
time  ot  a sloping  roof  is  adjusted  by  multiplying  the 
argument  of  the  hyperbolic  function  in  eq  8 by  the  cor- 
rection factor  F.  This  correction  is  included  in  the 
program  in  the  Appendix. 


The  two-dimensional  analysis  given  above  was  simpli- 
fied in  that  the  drain  was  assumed  to  be  a longitudinal 
gutter  capable  ol  being  an  infinite  sink.  Other  important 
factors  like  radial  (low,  drain  size,  basal  layer  permeability, 
and  gutter  slope  are  considered  below. 


RAINFALL  INTLNSITY-DUKATION  EFFECTS 

The  weight  added  to  a root  by  transient  water  in  both 
the  unsaturated  and  saturated  layers  increases  with  both 
the  intensity  and  duration  of  the  rainfall.  Snow-covered 
roofs  frequently  experience  rainfalls  which  would  cer- 
tainly cause  collapse  if  the  rainstorms  were  of  sufficient 
duration  to  allow  complete  soaking  of  the  unsaturated 
layer  and  lull  development  ot  the  saturated  layer. 

I ortunately,  the  duration  of  a rainstorm  generally  de- 
creases with  increasing  intensity  (Wisler  and  Brater  1965); 
hence  some  storm  of  intermediate  intensity  but  longer 
duration  may  actually  be  the  “design  basis  storm."  The 
fact  that  the  weight  of  the  unsaturated  layer  only  in- 
creases as  /l,,Jand  the  weight  of  the  saturated  layer  only- 
increases  as  I'12  suggests  that  both  duration  and  intensity 
are  important  parameters.  Accordingly,  we  might  sus- 
pect that  long  storms  of  less  intensity  may  add  more 
weight  to  some  roofs  whereas  others  may  be  loaded 
more  by  short,  intense  rainstorms. 

The  25-year  rainstorm  for  Hanover,  New  Hampshire, 
is  used  here  as  a design  basis  because  we  are  concerned 
only  with  the  coincidence  of  a heavy  rainfall  on  an  al- 
ready heavy  snow  load.  The  probability  of  a 25-year 
rainstorm  falling  on  a heavy  snowload  is  not  known 
but  this  example  does  provide  an  illustration.  If  it  is 
decided  that  a larger  safety  margin  is  needed,  a larger 
return  period  can  be  used,  but  because  the  roof  weight 
is  not  too  sensitive  to  the  storm  intensity,  a large  in- 
crease in  return  period  would  give  a small  increase  in 
roof  weight.  For  example,  increasing  the  return  period 
to  100  years  in  Hanover  would  only  increase  the  weight 
of  the  saturated  layer  by  10%.  However,  if  the  maxi- 
mum probable  precipitation  were  used,  the  weight  of 
the  saturated  layer  would  increase  by  about  80%  above 
that  calculated  for  the  25-year  return  period. 

From  Niedringhaus  (1973),  a 2.09  xl  O'5 -m/s  (3.0- 
in./h)  rainfall  lasting  1800s  can  be  expected  once  every 
25  years  in  Hanover,  N.l  I.  T rom  the  exponential 
formula  suggested  by  Wisler  and  Brater  (1965),  the 
intensity  / and  frequency  t of  25-ycar  rainfalls  in 
Hanover  arc  found  to  be 

i=  5.465  (r  + 360)-0-725  (mm/s)  (21) 


Duration  (s) 


Figure  9.  Total  live  load  ft',  fraction  of  rainfall  retained,  and  depth  of  water  impounded 
on  a flat,  two-dimensional  roof  for  the  25-year  rainfall  at  Hanover,  N.H.  The  maximum 
weight  occurs  for  a rainfall  of  about  75,000  s duration  hut  is  not  very  sensitive  to  dura- 
tion. 


\ 
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where  the  duration  r is  expressed  in  seconds.  The  total 
precipitation  i t increases  with  duration  so  that  lower 
intensity  rainstorms  always  deposit  more  water  on  the 
roof  than  the  higher  intensity  storms.  Clearly  some 
storm  of  intermediate  intensity  will  prove  to  be  the 
design  basis  rainstorm,  since  more  of  the  water  can 
drain  from  the  roof  during  the  lower  intensity  storms 
but  lower  intensity  storms  produce  more  total  precipi- 
tation. 

We  make  the  assumption  that  the  intensity'  of  the 
rain  is  constant  during  the  storm,  although  the  worst 
case  could  arise  if  the  intensity  were  unevenly  dis- 
tributed during  the  storm.  Unfortunately,  statistics 
on  the  intensity  distribution  are  not  available  and  no 
definitive  statements  about  the  effects  of  an  uneven 
distribution  can  be  made  from  the  equations  presented 
here.  If  there  is  a serious  concern  about  distribution, 
it  could  best  be  accommodated  by  increasing  the  de- 
sign basis  to  50  or  100  years. 

The  intensity-duration  effect  is  illustrated  by 
calculating  the  maximum  liquid  retention  of  the  satu- 
rated and  unsaturated  layers  as  separate  systems  and 
then  the  retention  of  the  entire  snow  load.  To  evalu- 
ate the  maximum  weight  of  the  saturated  layer  alone 
due  to  a 25-year  rainfall,  w'e  set  (/0(0)  equal  to  zero 
and  the  infiltration  / in  cq  8 equal  to  the  rainfall  /. 
Then  the  maximum  depth  of  water  impounded  on  a 
flat  roof  is 

d0  - i'12  h\]l2  L tanh  — — — t.  (22) 
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Tor  typical  parameters  for  a snow  with  a highly  perme- 
able basal  layer,  Figure  9 shows  that  the  depth  of  water 
ponded  on  the  roof  is  at  a maximum  for  an  intermedi- 
ate value  of  rainfall  duration.  When  the  rainfall  in  this 
example  lasts  100,000  s (27.8  h)  at  an  intensity  of 
0.00129  mm/s  (0.183  in./h),  the  depth  of  the  ponded 
layer  reaches  1 26  mm  (5  in.).  Only  51  % of  the  water 
reaching  the  saturated  layer  in  this  example  is  retained 
at  the  end  of  the  period  of  infiltration.  This  liquid 
would  add  53.4  kg/m2  (1 1.1  lb/fl2)  of  weight  or  about 
30%  of  a normal  design  load.  The  depth  of  water 
would  be  significantly  increased  if  1 ) the  roof  were 
covered  by  fresh  snow,  2)  the  culvert  or  drains  were 
incapable  of  handling  the  water  flowing  to  them,  or  as 
shown  below,  3)  the  flow  to  the  drains  occurred  as 
radial  rather  than  parallel  flow. 

The  maximum  weight  of  the  liquid  in  the  unsaturated 
layer  alone  can  be  determined  for  the  25-ycar  rainfall 
by  combining  cq  21  with  the  weight  of  the  liquid  in  the 
unsaturated  layer  which  is  given  by 

IFu=hpw0|(//aA'u)'/.t(i  -Swi)  + Sj  . (23) 

The  design  basis  storm  which  gives  the  maximum  in- 
crease in  the  weight  of  the  unsaturated  layer  is  calcu- 
lated for  snow  which  is  wet  but  has  no  antecedent 
flow  from 

h<5e  =0.149  (otu)’/3  r0-777  . (24) 
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I ho  u ilu-vil  dui  alion  would  be  longet  it  the  snow  were 
dry  prior  to  the  onset  ol  i a in  (this  option  is  included 
in  the  computer  program)  or  shortei  il  antecedent 
llow  existed  in  the  snow.  I lowever , the  weight  ol  the 
unsatui aled  layer  in  this  example  would  not  be  greatly 
affected.  I c|uation  24  clearly  shows  that  the  critical 
duration  ol  the  rainstorm  (and  the  weight  ol  the  un- 
saturated layer)  increases  rapidly  with  the  depth  ol 
the  snow  but  is  much  less  sensitive  to  the  permeabil- 
ity o!  the  snow.  I or  the  case  ol  Hanover  with  a snow 
ol  0.5-m  depth,  0.5  eltcctive  porosity  and  lO'^-m- 
permeability,  the  25-year  rainstorm  that  increases  the 
weight  ol  the  unsaturated  layer  to  a maximum  has  a 
duration  ol  8000  s (2.2  h)  and  an  intensity  ol  7.83x10'* 
inm/s  (0.23  in./h).  This  storm  adds  62.6  kg/m’  (12.8 
lb ' It-)  of  transient  liquid  weight  to  the  unsaturated 
layer  or  about  35%  ol  the  usual  design  live  load.  It  is 
very  important  to  note  that  about  18.9  kg/m2  (3.87 
Ih/ft - ) of  additional  liquid  is  retained  as  the  capillary 
liquid.  The  weight  of  the  ice  would  be  about  21  I 
kg/m2  (43.2  I b / f t - ) . Clearly,  the  additional  weight  of 
the  liquid  must  be  considered. 

In  these  examples,  the  unsaturated  layer  can  carry 
about  the  same  load  as  the  saturated  layer;  this  has 
important  implications  for  roofs  ol  all  slopes.  If  these 
weights  were  simply  added,  the  total  effect  would  re- 
quire an  increase  ol  the  design  load  of  more  than  two- 
thirds  to  account  for  the  total  weight  of  liquid  water 
on  the  roof.  However,  a more  complete  calculation  ol 
the  combined  effects  of  these  two  layers  is  necessary. 

I Iris  is  illustrated  here  for  a flat  roof  in  Hanover. 

f or  a flat  root  the  total  weight  of  the  liquid  per 
unit  area  is  given  by 

U Apw0l(//«fcu>l/JO  -Swl)  + Swi| 

+ 0.5  ncl0  0pw  (25) 

where  S’  is  the  itreducible  water  salutation  and 


rainstorm  lot  Hanover,  we  find  using  the  computer  pro- 
gram in  the  Appendix  that  the  maximum  weight  of  liquid 
carried  by  the  Hat  rool  increases  with  rainfall  duration 
until  about  75,000  s (20.3  h).  As  shown  on  I igure9,ihe 
total  weight  due  to  the  liquid  is  not  sensitive  to  the 
duration  ol  the  rainstorm  as  long  as  the  duration  t is 
well  in  excess  of  the  time  tu  necessary  to  penetrate  the 
unsaturated  layer.  In  (act,  the  weight  ol  the  transient 
liquid  increases  to  about  79  kg/m2  (16.4  lb/ft2)  before 
the  action  of  drainage  becomes  more  important  than  the 
increase  ol  total  precipitation  with  duration.  For  a rain- 
storm of  75,000-s  duration,  about  66%  of  the  rainfall 
is  retained  on  a llal  rool  by  the  end  ol  the  rainstorm. 

I or  this  situation  the  total  liquid  adds  98  kg/m2 
(20.4  lb,  It2)  ol  liquid  lo  the  snow  on  the  roof  ol  which 
the  capillary  liquid  (held  by  capillary  forces)  is  19  kg, 
m2  and  the  transient  liquid  is  79  kg/rn2.  The  design 
live  loads  on  most  roofs  would  have  to  be  increased  by 
about- 50%  to  account  for  the  total  liquid  weight.  I he 
total  live  load  in  this  example  reaches  304.3  kg/m2 
(62.4  lb/lt2);  thus  it  is  easily  seen  why  Thorburn  and 
Schriever  (1959)  were  concerned  about  roof  failures 
resulting  from  rain  tailing  on  a heavy  snowload. 

Because  of  the  large  number  of  parameters  affecting 
the  retention  of  rainwater  on  a roof,  a wide  range  ol 
calculated  loads  might  be  possible.  I lowever,  it  has 
been  shown  that  the  effects  of  snow  permeability  and 
rainfall  intensity  are  reduced  by  the  square  or  cube  root 
dependency,  and  as  suggested  by  Figure  9,  the  weight 
of  transient  liquid  on  a Hat  roof  is  not  very  sensitive  to 
duration  of  the  25-year  rainstorm.  Accordingly,  a wide 
variety  of  snow  conditions  and  rainstorms  can  produce 
about  the  same  amount  of  transient  roof  weight. 

The  drain  spacing  will  affect  the  saturated  layer,  and 
for  this  particular  example,  the  saturated  layer  is  very 
sensitive  to  the  slope  of  the  roof.  Other  important 
effects,  especially  the  radial  movement  of  water  to 
drains  and  gutter  overflow,  arc  introduced  in  the  next 
section. 


0 for  t < t , 


(1/akjV-  I tanh 


(r-/u)  for  T >tu. 


tu,  the  time  delay  before  the  infiltrating  water  reaches 
the  roof,  is  given  by 

'u  Wc  («V2)’,/3  • (27) 

I or  0.5-m  depth,  0.54  porosity,  10-m  drain  spacing, 
I0'l0-m2  permeability  for  unsaturated  snow,  10'9-m2 
permeability  for  saturated  snow,  and  the  25-year 


MISCELLANEOUS  EFFECTS 
Radial  flow  to  drains 

On  roofs  with  internal  drains,  radial  rather  than 
parallel  flow  can  occur  as  water  moves  over  the  rool 
to  the  drains.  Accordingly,  the  previous  analysis  lor 
flow  along  a two-dimensional  roof  must  be  extended 
to  allow  calculation  of  the  larger  weight  due  to  im- 
pounded water  associated  with  radial  llow  to  a drain. 

Steady  flow  q through  the  saturated  laser  on  a 
flat  roof  {0  0)  is  given  by 
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figure  10.  Saturated  layer  of  snow  on  a Hat  roof  shown  tor  the 
cases  where  the  saturated  layer  does  (!)  and  does  not  (H)  cover  the 
drain.  Note  that  the  effective  radius  of  the  drain  (ruJ  has  been  in- 
creased by  providing  a collector  space  between  the  drain  pipe  and 
the  grate. 


q = '(*k%2nrd^.  (28) 

The  infiltrating  water  reaching  the  saturated  layer  in- 
creases the  flux  toward  the  drain  according  to 

2nd  - ^ = 2w0  ^ . (29) 

Upon  combining  eq  28  and  29  to  eliminate  q and 
integrating  from  the  drain  ru  to  any  position  r tor 
steady  flow: 

^2=^-^('2-n5)  + ^lnW'D)  (30) 

where  d L is  the  thickness  of  the  saturated  layer  at  the 
perimeter  of  the  drain,  rD  is  the  drain  radius,  and  re 
is  the  roof  radius  (see  Fig.  10).  The  size  and  spacing 
of  the  drains  are  critical,  since  the  drains  must  be 
capable  of  removing  all  the  water  without  it  ponding 
over  them.  As  shown  on  Figure  10,  the  weight  of  the 
saturated  layer  is  much  greater  when  the  drain  is  in- 
capable of  accepting  all  the  water  moving  through  the 
saturated  layer.  The  vertical  percolation  through  the 
saturated  layer  directly  over  the  drain  is  independent 
of  the  thickness  of  that  layer  because  the  water  pres- 
sures on  both  the  upper  and  lower  surfaces  arc  equal. 
Therefore  the  flow  through  the  drain  QD  is  given  by 
the  inequality: 

|0D|  • aA \nrl  . (31) 

The  total  flow  reaching  the  saturated  layer  is  nr  2 1. 

Thus 

(rjr0)2f-  cd\  (32) 


{rjru)2  I can  be  mote  than  an  order  of  magnitude 
greater  than  oA’s;  however,  as  noted  before,  the  time 
required  for  the  saturated  layer  to  develop  is  very  much 
greater  than  the  duration  of  intense  rainfalls.  Neverthe- 
less there  are  several  important  conclusions  to  be  drawn 
at  this  point. 

1 . The  rate  of  flow  to  the  drain  is  very  sensitive  to 
the  ratio  of  drain  spacing  to  draiit  size  as  shown  on  the 
left  hand  side  of  eq  32.  This  ratio  can  be  reduced  by 
increasing  the  size  of  drains  and/or  decreasing  the 
spacing  between  the  drains.  Since  internal  drains  are 
expensive  and  lend  to  be  troublesome,  it  is  important 
to  note  that  the  effective  diameter  of  a drain  can  be  in- 
creased by  simply  increasing  the  size  of  the  grating  over 
the  drain.  The  area  available  for  flow  to  the  drain  in- 
creases as  r 5 for  Case  I (see  Fig.  1 0)  and  as  rD  for 
Case  II.  By  increasing rD  the  increased  weight  due  to 
ponding  over  the  drain  can  be  greatly  reduced. 

2.  One  severe  problem  with  drains  is  the  possibility 
of  icing  because  of  freeze-thaw  cycles  and/or  the 
densification  of  the  saturated  snow  overlying  the  drain. 
When  the  snow  overlying  the  grating  is  wet  it  will  retain 
by  capillarity  approximately  20  mm  (0.8  in.)  of  liquid 
water  in  a saturated  layer  at  the  very  base  of  the  snow. 

If  weather  conditions  change  and  this  layer  refreczes, 

it  becomes  a nearly  impermeable  "ice  layer."  If  a subse- 
quent rainstorm  occurs,  discharge  into  the  drain  is 
impeded  by  the  reduced  permeability  of  this  layer  and 
more  of  the  rainwater  could  be  retained  on  the  roof. 

The  “ice  layer”  undergoes  a thermodynamic  meta- 
morphism  in  the  presence  of  liquid  water  which  eventu- 
ally restores  the  permeability  of  the  snow  over  the 
drain;  however,  it  takes  several  days  to  cause  a signifi- 
cant increase  in  As  a result  the  weight  added  to  a 
roof  by  a steady  rainfall  lasting  several  days  could 
approach  the  total  weight  of  the  rainwater  falling  on 
the  roof.  One  obvious  approach  to  solving  this  problem 
is  to  heal  the  grating  in  order  to  melt  away  the  20-mm 
ice  layer.  1 his  solution  can  be  partly  accomplished  by 


and  there  is  no  ponding  over  the  drain  as  long  as 
(rc,V| j)2  / is  less  than  aA’s.  It  is  conceivable  that 
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figure  II.  rhree  advantages  of  a recessed  drain 
over  a flush  drain:  / ) the  layer  of  snow  saturated 
by  capillary  retention  will  he  limited  to  the  drain 
itself,  thus  reducing  the  root  weight,  2)  this 
saturated,  capillary  layer  will  be  less  likely  to 
free/e  and  block  the  drain,  and  3)  the  lowered 
drain  creates  a "tension  gradient" in  the  snow 
which  slightly  increases  the  rate  of  flow  to  the 
drain. 

recessing  the  grating  as  shown  on  Figure  1 1 . A recess 
oi  50  mm  (1.97  in.)  has  several  advantages  over  the 
flush  grating  shown  in  Figure  10.  First,  the  layer  of 
snow  saturated  by  capillary  retention  is  limited  to 
the  drain,  thus  reducing  the  weight  that  this  saturated 
layer  exerts  on  the  entire  roof.  Second,  internal 
drains  form  a natural  heat  leak  from  the  building 
which  should  help  to  prevent  the  recessed  snow  from 
freezing.  The  greater  the  heat  leak,  the  less  the  chances 
for  ice  layer  formation  to  cause  drain  blockage. 

3.  The  pressure  in  the  water  at  the  lowest  point  in 
the  snow  is  approximately  atmospheric  and  decreases 
with  height  above  that  point.  This  establishes  a 
“tension  gradient"  which  has  been  ignored  in  the 
derivation  of  the  flow  equations  given  here  but  in  fact 
would  increase  the  flow  rate  toward  a recessed  drain. 

4.  Assuming  the  drain  is  not  blocked  and  is  capable 
ol  handling  the  infiltrating  water  without  forming  a 
ponded  layer  directly  over  the  drain,  the  profile  of 
water  in  the  saturated  layer  is  given  by 

(^p)'  - “ (2  In  R/R0  -ff2 +/?£)>*  (33) 

' ' rc 

where  re  is  the  radius  of  the  roof, 

«o  rD  /'e  <34> 

and 

R r/rc  . (35) 

The  dimensionless  ratio  (2cdrJf)'/!  d/r e is  shown  as  a 
function  of  R for  various  values  of  in  Tigurc  12. 
The  effect  of  drain  size  is  shown  by  the  large  increase 


in  the  thickness  ol  the  saturated  layer  for  small  values 
o\Ru. 

5.  On  a roof  with  radial  flow,  a larger  volume  of 
water  is  retained  as  a saturated  layer  than  on  a roof  with 
parallel  flow  to  a drain.  To  illustrate  this  point  the  pro- 
file of  the  saturated  layer  on  a two-dimensional  toot  is 
given  by 

(2^r)  f =\/2  (1  -Rl)Vl  (36) 

where 

Rt=x/L  (37) 

This  is  also  shown  in  Figure  1 2 for  a drain  size  of  0.1 
of  the  roof  length.  The  volume  of  ponded  water  on  a 
root  with  radial  flow  is  given  by 

r,, 

V = 2tt0  / rddr.  (38) 

J0 

Hence  the  ratio  of  the  volume  of  ponded  water  for 
radial  flow  to  that  for  parallel  flow  can  be  calculated 
as  a function  ol  R^ . This  ratio  is  shown  as  a function 
of  R u in  Figure  13  for  roofs  of  equal  area.  On  the 
radial  roof  in  this  figure,  there  is  a single,  centrally 
located  drain,  and  the  two-dimensional  roof  has  a f ixed 
drain  size.  For  R [(  equal  to  0.1 , the  drain  is  '/,0  the 
size  of  the  roof  and  the  saturated  layer  for  radial  flow 
has  78%  more  weight  than  for  parallel  flow . For  R 
equal  to  0.01  which  is  probably  a more  common  value, 
the  saturated  layer  for  radial  flow  contains  1 86'  more 
water  than  for  parallel  flow.  The  advantages  of  in- 
creasing the  relative  size  of  the  drain  R0  lor  radial 
flow  by  using  an  enlarged  grating  over  the  drain  as 
suggested  earlier  have  now  been  verified.  Likewise, 
the  advantage  of  draining  a roof  by  parallel  flow  to 
a gutter  rather  than  radial  flow  to  drains  has  been  es- 
tablished on  a quantitative  basis. 

As  for  parallel  flow,  the  rate  of  increase  ol  the 
saturated  layer  in  radial  flow  must  be  considered.  I he 
thickness  of  the  saturated  layer  at  the  perimeter  of  the 
roof  increases  for  a constani  value  of  infiltration  accord- 
ing to 

d0  = d0( 0)  + zl  (l/ak)  tanh  | (a kl)Vl  t/pA  | (39) 

where 

zl  = \ /r\  In  re/ru  - \ : r*  -*  r£  . (40) 


Figure  12.  Fro  file  ol  the  saturated  layer  (d/rc)  shown  as  a function  of  R for  various  drain 
sizes  (R0J  for  radial  How  over  a flat  roof.  The  effect  of  decreasing  the  drain  size  is  a rapid 
increase  in  the  thickness  ol  the  saturated  layer.  For  comparison,  the  corresponding  profile 
is  given  for  the  two-dimensional  roof  considered  previously  (dashed  line).  The  thickness 
of  the  saturated  layer  is  somewhat  greater  for  radial  than  for  parallel  flow. 
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figure  13.  Ratio  of  the  volume  of  water  in  the  saturated  layers  of  radial  and  two- 
dimensional  roofs  shown  as  a function  of  the  drain  size  on  a radial  roof. 
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Two  d'mens»onoi  Roof 


Snow 


a.  Two-dimensional  roof  with  a recessed  gutter  carries  the  minimum 
weight  when  the  gutter  depth  exceeds  the  water  depth. 


b.  When  the  gutter  cannot  move  all  of  the  water  to  the  drain  without 
overflowing,  the  roof  slope  is  effectively  decreased  and  the  weight  due 
to  impounded  water  is  greatly  increased. 


Figure  14.  Flow  on  a flat  roof  for  two  conditions. 


From  eq  38  the  maximum  possible  volume  on  the  roof 
is 

V = V2  r,g>  \/l/ak  r\ 

\/|n  R2  -R2  + R^~  In  Rfr  RdR  . (41) 

''D 

In  the  formulation  of  the  computer  program  it  is 
assumed  that  the  volume  of  the  saturated  layer  in- 
creases directly  as  dQ  (t). 

Flow  along  gutters  on  snow-covered  roofs 

We  have  shown  that  planar  two-dimension  ll  roofs 
tend  to  carry  less  weight  in  the  saturated  layer  than 
radial  roofs.  The  advantage  of  two-dimensional  roofs 
is  greatly  reduced  when  the  longitudinal  gutter,  which 
collects  and  moves  the  water  to  the  drains,  cannot 
move  the  water  to  the  drains  without  overflowing. 

That  is,  when  the  thickness  of  the  layer  of  water  in  a 
recessed  gutter  exceeds  the  depth  of  the  gutter,  the 
overflowing  gutter  effectively  reduces  the  slope  of  the 
roof  in  order  to  increase  its  own  slope.  As  shown  in 
Figure  14,  this  phenomenon  causes  a large  increase  in 
the  weight  on  the  roof  because  small  decreases  in  roof 
slope  produce  large  increases  in  the  volume  of  the 
ponded  water  (sec  Fig.  5).  For  level  gutters,  the  depth 


of  the  saturated  layer  in  the  gutters  assumes  the  elliptical 
profile  of 


Flcnce  the  overflow  would  most  likely  occur  near  the 
roof  edge  and  extend  along  most  of  the  length  of  the 
gutter.  For  a steeply  dipping  gutter,  the  depth  of  the 
saturated  layer  increases  linearly  toward  a drain  and 
reaches  a maximum  just  short  of  the  drain.  In  this 
case  the  additional  weight  is  concentrated  near  the 
drain  so  that  an  overflow  would  nut  cover  such  a 
large  area  of  the  roof.  Gutters  with  gradual  slopes, 
like  roofs  with  gradual  slopes,  lie  somewhere  between 
the  two  extremes. 

Obviously  much  can  be  gained  by  increasing  the 
slope  of  a gutter.  Unfortunately  it  is  difficult  to  solve 
the  governing  equations,  and  unlike  the  case  for  a 
shallow  roof  where  the  water  volume  was  used  as  the 
dependent  variable  instead  of  depth,  we  must  deal 
with  depth  explicitly  since  the  depth  of  the  water  in 
the  gutters  is  the  limiting  criterion  for  design.  The 
volume  of  water  in  the  gutters  per  sc  contributes  very 
little  to  the  roof  weight.  Ideally  the  gutters  should  be 
able  to  handle  all  the  rain  intercepted  by  the  rool 
without  overflowing  but  this  is  very  unlikely  for  nor- 
mal roof  design.  Sloping  the  gutters  reduces  the  weight 
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addition  by  overflow  somewhat,  but  even  so,  the 
gutters  would  have  to  be  extremely  large  to  handle 
all  the  drainage,  I or  a Hat  culvert,  the  maximum 
depth  ol  water  r/0  in  the  gutter  is  given  by 


(?)W 


where  is  the  length  ol  the  gutter,  21.  is  the  width 
ol  the  rool  and  u is  the  width  ol  the  gutter.  Equation 
42  shows  that  the  depth  ol  water  in  the  gutter  in- 
creases with  the  spacing  between  the  drains  and  as 
the  square  root  ol  the  ratio  ol  rool  to  gutter  widths 
(2 L/w)'  . For  reasonable  values  of  (2L/w)' ' and 
(I/a A’s):-',  the  depth  of  the  gutter  must  be  excessively 
large  unless  the  drain  spacing  is  no  more  than  a meter. 
Accordingly  a recessed  gutter  similar  to  the  suggested 
drain  design  (see  Fig.  1 1 ) would  be  advantageous.  The 
main  advantage  ol  this  design  is  that  the  longitudinal 
flow  along  the  gutter  occurs  in  the  snow-free  space 
beneath  the  grating;  hence,  flow  along  the  gutter  would 
occur  much  more  readily.  For  given  conditions  the 
tlow  rate  would  be  greatly  increased  at  small  values  of 
water  depth.  Accordingly  overflow  would  not  in- 
crease the  roof  weight  by  inhibiting  llow  to  the  gutter. 

Snow  structure 

In  most  areas  of  heavy  snowfall,  ice  layers  tend  to 
form  in  the  snowcover  because  of  the  preferential  re- 
tention of  liquid  water  at  certain  horizons.  These  ice 
layers  interrupt  the  downward  movement  of  water  in 
the  unsaturated  zone  causing  ponded  layers  of  water 
to  form  above  them.  In  extreme  cases  these  ice  layers 
almost  completely  eliminate  downward  water  move- 
ment and  it  would  be  necessary  to  manually  penetrate 
such  an  ice  layer  in  the  event  of  heavy  rains. 

Generally  ice  layers  arc  discontinuous  and  do  not 
provide  a complete  interruption  of  the  vertical  seepage. 
In  effect,  discontinuous  ice  layers  cause  a stepwise  flow 
of  water  down  to  the  roof.  There  are  two  important 
effects  of  this  mode  of  flow.  First,  more  than  one 
saturated  horizon  will  exist  since  each  ice  layer  as  well 
as  the  roof  tend  to  cause  v/atcr  ponding.  Second, 
some  areas  of  the  snowcover  may  be  bypassed  and  re- 
main relatively  dry  while  the  flow  is  concentrated  in 
zones  of  high  water  saturation  and  high  flow  rates. 
Accordingly,  we  expect  an  uneven  distribution  of  the 
weight  on  the  roof  due  to  the  liquid  water. 

It  is  difficult  to  generalize  about  the  relative  roof 
loads  for  homogeneous  snow  vs  snow  containing  ice 
layers.  This  problem  is  compounded  by  the  highly 
variable  nature  of  ice  layers  in  time  and  space  and  by 
the  lack  of  studies  of  the  particular  properties  of  snow- 


on  roots.  1 he  extent  and  nature  ol  ice  layers  in  the 
snowcover  on  the  ground  change  rapidly  during  the 
warmer  spring  months  when  heavy  snowloads  associ- 
ated with  heavy  rains  are  most  likely  to  occur.  In 
general  it  seems  sale  to  assume  that  the  presence  ol  ice 
layers  will  not  significantly  increase  the  weight  due  to 
the  presence  ol  transient  water.  However,  it  must  be 
emphasized  that  experimental  observations  of  the  de- 
tailed structure  of  snowcovers  on  roots  are  necessary 
to  test  the  validity  of  generalizing  from  snow  on  the 
ground  to  snow  on  roofs. 

Basal  layer 

There  are  two  characteristics  of  the  basal  layer  of 
snow  which  are  important  the  basal  layer  permea- 
bility and  open  channels.  It  was  shown  earlier  that  the 
permeability  of  the  basal  layer  is  ol  reduced  importance 
because  it  enters  the  equations  as  a square  root.  I urthcr- 
more',  at  small  values  of  time 


ML  t * 
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Therefore,  from  eq  16  for  small  values  of  time 

do  =■  d0(0)  + m (44) 

which  shows  that  initially  the  thickness  of  the  saturated 
layer  increases  linearly  with  time  and  is  relatively  insensi- 
tive to  the  permeability.  This  is  illustrated  on  ( igurc  1 5 
for  two  values  of  permeability.  Note  that  the  effect  of 
permeability  increases  with  time. 

The  higher  permeability  used  in  this  example  ( 1 O' 1 
mm2)  could  only  represent  a well-soaked  snow  with 
some  open  channels  formed  by  melting  in  this  basal 
layer.  The  maximum  thickness  of  the  saturated  layer 
is  10  times  larger  for  the  snow  of  lower  permeability 
and  its  time  response  is  much  slower.  These  are  both 
important  considerations  regarding  the  weight  of  the 
basal  layer.  More  observations  of  roof  snow  are  needed 
to  test  the  permeability  values  assumed  here. 


CONCLUSION 

Heavy  rains  tailing  on  roofs  carrying  a heavy  snow 
load  are  an  important  aspect  of  roof  failures.  The 
weight  of  the  water,  both  capillary  and  transient  liquid, 
should  be  considered  as  an  additional  load  which  must 
be  supported  by  the  roof.  I he  maximum  weight  ol 
liquid  depends  on  properties  ol  the  snow,  on  size,  slope 
and  shape  of  the  roof,  on  spacing  and  type  ol  the 
drains,  and  on  intensity  and  duration  ol  the  rainstorms. 
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Figure  1 5.  Depth  of  the  saturated  layer  (dj  shown  as  a function  of  time  for  an 
intense  rainstorm  for  two  values  of  permeability  of  the  basal  layer.  The  maximum 
depth  obtained  differs  by  a factor  of  10,  although  the  permeabilities  differ  by  a 
factor  of  / 00.  The  more  permeable  basal  layer  has  a much  more  rapid  time  re- 
sponse. 


I 
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While  the  depth  of  the  snow  is  very  important,  the 
effect  of  permeability  of  the  snow  is  reduced  by  its 
square  or  cube  root  dependence.  However,  the 
presence  of  open  channels  at  the  base  of  the  snow 
will  have  an  effect  on  reducing  the  weight  carried  by 
the  roof.  1 he  roof  size,  especially  for  roofs  with 
internal  drains,  has  a very  large  effect,  since  the  pond- 
ing of  water  over  the  rpof  increases  as  the  square  of 
the  drain  spacing.  The  effect  of  giving  the  roof  even 
a small  inclination  is  very  large  and  must  be  consid- 
ered in  order  to  reduce  the  depth  of  the  ponded  layer. 
Lnlarging  and  recessing  the  drains  increase  the  rate  of 
drainage  and  reduce  the  chances  of  blockage  by  re- 
freezing. 

Lven  for  steeply  dipping  roofs,  much  water  can  be 
retained  in  the  unsaturated  snow.  This  water,  both 
capillary  and  transient,  must  be  considered  since  one 
example  for  a 25-ycar  rainstorm  in  Hanover,  New 
I lampshirc,  shows  that  61  kg/m2  (12.5  lb/ft2 ) of 
weight  could  be  added  to  the  roof  by  the  liquid  in  the 
unsaturated  snow  alone.  The  25-year  rainstorm  was 
chosen  for  use  a-,  an  example  in  order  to  illustrate 
the  effect  of  intensity  vs  duration.  If  intense  rain- 
storms lasted  long  enough  they  would  certainly  cause 
many  failures  of  flat  or  shallow  roofs  covered  by  snow, 
fortunately,  more  intense  storms  terminate  before  the 
ponded  layer  can  be  fully  developed  and  less  intense 
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storms  cause  only  a shallow  ponded  layer  to  develop 
because  of  concurrent  drainage. 

For  the  example  given  lor  Hanover  the  25-year  rain- 
storm falling  on  a heavy  snow  load  could  increase  the 
total  roof  load  due  to  liquid  water  to  98  kg/m2  (20 
lb/ft2).  The  normal  design  load  would  have  to  be  in- 
creased by  about  50%  to  account  for  this  effect.  If  the 
snow  load  is  unevenly  distributed  (see  cover  photograph), 
the  local  load  due  to  rain  falling  on  the  snowdrift  can  be 
as  large  as  that  produced  by  covering  the  entire  roof  with 
the  drift. 
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APPENDIX  A.  COMPUTER  PROGRAM  CALCULATING 
ROOF  LOADS  FROM  RAIN-ON-SNOW 


\ 

i 


This  computer  program  in  Dartmouth  Basic  (see 
Pig.  Al ) calculates  total  weight,  traction  ol  rainfall  re- 
tained, weight  ol  the  saturated  layer  and  weight  of  the 
unsaturated  layer  as  a function  of  duration  lor  a design 
basis  rainstorm  falling  on  a snow-covered  roof.  I he 
input  required  to  run  the  program  is  as  follows  (in 
MRS  units): 

1 ) Specify  A,  ti,  C 

These  are  the  constants  in  a design  basis  rainstorm 
described  by 

/ - /I  (r  + B)c  (Al) 

and  must  be  obtained  from  climatological  records  for 
each  area. 

2)  Specify  //,  L,P,S,K],  K2 

H is  the  thickness  ol  the  snow.  L is  the  distance  be- 
tween the  parapet  and  drain  for  parallel  flow  to  gutters 
or  the  roof  edge,  or  the  shortest  distance  between  the 
parapet  and  drain  for  radial  flow.  / 1 is  the  porosity  of 
the  snow.  5 is  the  irreducible  water  saturation  of  the 
snow.  K\  and  K2  are  the  intrinsic  permeabilities  of 
the  saturated  basal  layer  and  the  unsaturated  upper 
layers,  respectively. 

3)  Specify  radial  (R)  or  parallel  [R)  flow 

Radial  flow  is  flow  to  a central  drain  and  parallel 

llow  is  llow  to  a gutter  which  carries  the  water  to 
drains.  It  is  assumed  that  a square  area  is  drained  in 
radial  flow  and  that  the  gutter  extends  the  entire  length 
of  the  roof  for  parallel  flow.  Other  configurations 
would  require  some  modification  of  the  program  it 
they  could  not  be  represented  by  these  assumptions. 

4)  Specify  drain  radius  (meters) 

The  effective  radius  of  the  drain  for  radial  flow. 

No  drain  or  gutter  overflow  is  assumed. 

5)  Specify  temperature  (”C) 

The  average  snow  temperature  is  required.  II  the 
snow  is  already  at  O 'C.,  it  is  assumed  that  only  the 
capillary  water  is  present  in  the  snow  such  that  no 
antecedent  flow  occurs  prior  to  the  onset  ot  the  rain- 
storm. 


6)  Specify  roof  slope  (degrees) 

The  roof  slope  is  input  in  degrees  of  slope. 

7)  Specify  start  lime,  stop  time,  time  step 

The  program  makes  calculations  for  each  value  of 
rainstorm  duration  chosen  here. 

8)  Another  time  step? 

The  operator  can  optimize  the  time  step  to  minimize 
computer  time  by  using  a small  time  step  at  lower 
values  ol  duration  and  increasing  the  time  step  at 
larger  values  of  duration  where  the  output  changes  more 
slowly. 

The  output  from  the  program  is  the  average  weight 
per  unit  area,  the  maximum  weight  per  unit  area 
occurring  anywhere  on  the  roof,  the  fraction  of  the 
total  rainfall  retained,  the  average  weight  of  the  satu- 
rated layer,  and  the  average  weight  of  the  unsaturated 
layer  for  each  specified  value  ol  rainstorm  duration. 
Other  outputs  such  as  depth  of  the  saturated  layer, 
weight  ol  the  capillary  water,  weight  of  refrozen  water 
or  fraction  of  wetting  front  penetration  could  be  made 
with  minor  modifications.  II  the  drains  or  gutters  were 
not  capable  of  removing  the  water  without  overflowing, 
some  modifications  would  have  to  be  made  to  the  pro- 
gram to  account  tor  the  additional  weight. 


LIST 


RFDTROS  02  MAR  77  13: 03 

10  REM  ROOF  WEIGHT  DUE  TO  RAIN  ON  SNOW*  MKS  UNITS 
100  PRINT  "SPECIFY  AfBfC" 

200  INPUT  AfBfC 

300  PRINT  "SPECIFY  HfLfPfSfK1fK2" 

310  INPUT  HfLfPfSfK1fN2 
TOO  LET  W1=917*H#< 1-P) 

405  PRINT  "RADIAL (R)  OR  PARALLEL <P>  FLOW" 

410  INPUT  W* 

415  IF  W*="P"  THEN  500 

417  IF  W$= "R " THEN  425 

418  GO  TO  405 

425  PRINT  "DRAIN  RADI US ( METERS ) " 

430  INPUT  G 

435  LET  J=.887*G/L 

440  LET  M= J/100 

445  LET  J1 = J+M 

475  FOR  R- J1  TO  1 STEP  M 

480  LET  N1=<2*L0G<R)-R~2+J~2-2*L0G< J) ) ~ < ,5>*R*J 
485  LET  N=N+N1 
490  NEXT  R 

500  PRINT  "SPECIFY  TEMP  < DEG  C>* 

510  INPUT  T9 

515  LET  W5=5. 73*T9*H*<P-1 ) 

525  PRINT  "SPECIFY  ROOF  SLOPE ( DEG ) " 

530  INPUT  Q 

540  LET  F= ( l+Q/2 .2)~(1.06) 

570  PRINT  "SPECIFY  START  TIME  f STOP  TIME  f TIME  STEP" 

580  INPUT  T2  f T3  f T4 

590  PRINT  "DURATION" f "TOT  W*f  "RETENTION" f "W(SAT)"f  "W(UNSAT)" 

600  FOR  T=T2  TO  T3  STEP  T4 
700  LET  I=A*(T+B)"C 
770  IF  T9<0  THEN  815 
780  IF  T9-0  THEN  800 
790  GO  TO  500 

800  LET  T1=<H*P*( 1-S) )/( < 5470000*K2# < I"2> ) " < 1/3) ) 

805  GO  TO  900 

815  LET  T1=H*P*((I/(5470000*K2)'"(1/3)*(1-S)+S)/I  + 0.001*U5/I 
900  LET  U=< <5470000*K1#I )~< 1/2) )*<T-T1 )/(L*P) 

910  LET  U^U*F 

1000  LET  D-0 

1005  IF  T<T 1 THEN  1 326 

1010  IF  W*="P"  THEN  1200 

1060  LET  N2=15.95*P#L*N*(I/<5470000*K1> >”< .5) 

1070  LET  Ul  = < <5470000*K1*I  ),'.5*(T-T1  ) ) / ( P*  ( L"2*L0G  < L/G ) - . 5*L~2+  . 5*G"2 ) " . 5 

1071  LET  U1-~F*U1 

1080  LET  W2~N2#(EXP(U1 )-EXP(-Ul ) )/<EXP<Ul )+EXP<-Ul ) ) 

1081  LET  W2=U2/<f1*<5470000*K1/I>~(.5>*<F-1>+1> 

1085  LET  D~.00101#W2/P 

1090  GO  TO  1322 

1200  LET  D-( ( I/(5470000*K1 ) )"< 1/2) ) *L* < ( EXP < U ) -EXP ( -U ) ) / ( EXP ( U ) +EXP ( -U ) ) ) 
1320  LET  W2"785*P*D/((F-1)*.1#(5470000*K1/I)'~(.5)+1) 

1322  LET  E~1 

1324  IF  T>T 1 THEN  1340 

1326  LET  W 2 0 


Figure  A 1.  Computer  program. 
18 


.?  • 


r 


1330  LET  E-T/ri 

1340  LET  W3» 1000*H*P*< ( . 005676# < I/K2 ) " . 3333 ) * ( 1~S ) ) 
1350  LEI  W 4 ™ 1 0 00 # H * P * S 
1355  IF  IV  0 THEN  1370 

1360  I El  W- -W.I  +W2+E*<  1-D/<H*F>  >*<W3+W4) 

1364  LEI  0 ■:  ( 1-D/<H*F>  >*E*W3+W2>/<  1000*1*7  ) 

1366  GO  10  1495 

1370  LET  W Wl+W2+E*< 1-U/(H*F) ) * < W3+W4 > +E4W5 

1380  LET  M • < LI.:  +1  *<  1-D/<H*F>  ) # < W3+W4 ) +E*W5 )/ < 1000* I #T> 

1495  LET  Y-  < 1-H/(H#F)  )*E*W3 

1 499  IF  T7:  0 THEN  1750 

1500  PRINT  TfW»V?W2»Y 
1600  NEXI  T 

1650  PRINT  " ANOTHER  TIME  STEP?* 

1660  INPUT  Z% 

16  70  IF  2*-.: -YES*  THEN  570 
1680  IF  ZI  "NO''  THEN  1700 
1690  GO  TO  1650 
1700  IF  E:  1 THEN  2000 

1710  PRINT  -DURATION  OF  MAXIMUM  TOTAL  WEIGHT?* 

1720  INPUT  T7 

1730  LET  T~ T7 

1740  GO  TO  700 

175 0 L E T W 6 ■ W 1 + 1 0 0 04  D* P + Y + W 4 

1800  IF  190  THEN  1950 

1900  LET  W6-W6+E4W5 

1950  PRINT  "MAXIMUM  TOTAL  WEIGHT  PER  UNIT  AREA  =*»W6 
l 2000  END 

R E A D Y 
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